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Abstract 

We describe a new technique which enables one to construct an SO(10)-symmetric 
fermion mass matrix model at the supersymmetric grand unification scale directly from 
the fermion mass and mixing data at the low energy scale. Applications to two different 
neutrino mass and mixing scenarios are given. 
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1 Description of the Technique 

A new bottom-to-up approach proposed by us is summarized briefly in this talk with 
more details and references presented elsewhere. We shall restrict our attention here to 
the construction of complex symmetric mass matrices arising with Higgs in the 10 and 126 
representations of SO (10). The procedure allows us to construct mass matrices which exhibit 
as simple an SO(10) structure as possible with the maximum number of texture zeros allowed 
for the neutrino mass and mixing scenario in question. 

1) Start from a set of quark and lepton masses and mixing matrices completely specified 
at the low scales. 

2) Evolve the masses and mixing matrices to the GUT scale by making use of the one-loop 
renormalization group equations (RGEs) for the minimal supersymmetric standard 
model. 

For this purpose we set Asusy = 170 GeV and A GUT = 1.2 x 10 16 GeV. Following Naculichjg 
we use the approximation that only the top and bottom quarks as well as the tau lepton 
contribute to the non-linear Yukawa terms in the RGEs. With a physical top mass expected 
to be near 160 GeV, we take the running mass to be m t (m t ) = 150 GeV and adjust m&(m;,) 
and tan/3, so consistency is achieved at Aqut which requires complete Yukawa unification 
with tan/3 ~ 48.9. We are working under the assumption that only one SO(10) 10 of Higgs 
contributes to the 33 elements of the up, down and charged lepton mass matrices. 

3) Construct a numerical set of complex-symmetric mass matrices M u , M D , M E and 
M N *ft = M Nmrac Mft 1 M NmracT for the up and down quarks, charged leptons and light 
neutrinos by making use of a procedure due to Kusenko, || now applied to both quarks 
and leptons. 
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Since the quark CKM mixing matrix is unitary and represents an element of the unitary 
group U(3), one can express it in terms of one Hermitian generator of the corresponding 
U(3) Lie algebra times a phase parameter a by writing 

V C km = U' L Ul = exp(iaH) (la) 

where 

iaH = ^ log v k ) — ^ — 16 

in terms of the eigenvalues Vj of Vckm, by making use of Sylvester's theorem. The trans- 
formation matrices from the weak to the mass bases are given in terms of the same generator 
but modified phase parameters such that 

U' L = exp(iaHx q ), Ul = exp [iaH(x q — 1)] (lc) 

and relation (la) is preserved. The complex symmetric quark mass matrices in the weak 
basis are then related to those in the diagonal mass basis by 

M u = 1J2d u u2 T } M d = UlD D Ul T (Id) 

It suffices to expand Vckm, U' l and Ul to third order in a in order to obtain accurate 
numerical results for the mass matrices M u and M D . Similar expressions can be obtained 
for the light neutrino and charged lepton mass matrices from the lepton mixing matrix and 
its eigenvalues with xt replacing x q . 

The parameters x q and xe control the choice of bases for the quark and lepton mass 
matrices, respectively. The up quark mass matrix is diagonal for x q = 0, while the down 
quark mass matrix is diagonal for x q — 1; likewise, the light neutrino mass matrix is diagonal 
for xi = 0, while the charged lepton mass matrix is diagonal for xe = 1. 

4) Vary x q and xg and the signs of the mass eigenvalues to search for simplicity in the 
SO (10) framework, i.e., pure 10 or pure 126 contributions for as many mass matrix 
elements as possible. 



For this purpose we note that 10's contribute equally to the down quark and charged lepton 
mass matrices, while 126 contributions differ by a factor of -3; likewise for the up quark and 
Dirac neutrino mass matrices. In terms of the Yukawa couplings and the appropriate VEVs, 
the mass matrices are given by 

M u =j:j^) Vm + E 3 f {12 ^w UJ 
M N Dirac =j2j(^) Vui -3j: j f^w uj 

M E = Ei/ (10^ d ._ 3E . / ( 1 26,) Wd . 

5) For the best choice of x q and xt which maximizes the simplicity, construct a simple 
model of the mass matrices with as many texture zeros as possible. 

6) Evolve the mass eigenvalues and mixing matrices determined from the model at the 
SUSY GUT scale to the low scales and compare the results with the starting input 
data. 

2 Application to Two Different Neutrino Mass and Mixing Scenarios 

We now illustrate the technique by applying it to two different neutrino scenarios, 
both of which explain the solar neutrino depletion data with the non-adiabatic Mikheyev- 
bmirnov- Wolfenstein@ (MSW) effect, where one includes the atmospheric neutrino depletion 
phenomenon U while the other accepts the cocktail model interpretation of missing dark 
matter. 

We take the same quark input data for both models. For the light quark masses, we 
shall adopt the values quoted by Gasser and Leutwyler|7| while the heavy quark masses are 
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specified at their running mass scales: 

m u (lGeV) = 5.1 MeV, m d (lGeV) = 8.9 MeV 

m c (m c ) = 1.27 GeV, m s (lGeV) = 175 MeV (3a) 

m t (m t ) = 150 GeV, m b (m b ) ~ 4.25 GeV 

For the Cabibbo-Kobayashi-Maskawa (CKM) mixing matrix, we adopt the following central 
values at the weak scale 

/ 0.9753 0.2210 (-0.283 - 0.126i) x 10~ 2 \ 

Vckm = -0.2206 0.9744 0.0430 (36) 

\0.0112 - 0.0012i -0.0412 - 0.0003i 0.9991 J 

where we have assumed a value of 0.043 for V cb and applied strict unitarity to determine 
V ub , V td and V ts . 



2.1 Lepton Scenario (A) involving the Non-adiabatic MSW Solar and Atmo- 
spheric Neutrino Depletion Effects 



In this scenario, we single out the central points in the two neutrino mixing planes 



8m\ 2 



5 x 10~ 6 eV 2 , 



sin 2 2^2 = 8 x 10~ 3 
2xlO~ 2 eV 2 , sin 2 2^3 =0.5 

With the neutrino masses assumed to be non-degenerate, we take for the lepton input 



om 2 3 



(4) 



and 



0.5 x 10~ b eV, 
0.224 x 10^ 2 eV, 



m 
m 



77V = °- 141 eV > 



m 1 



0.511 MeV 
105.3 MeV 
1.777 GeV 



LEPT 



0.9990 0.0447 (-0.690 - 0.310?) x 10~ 2 \ 

-0.0381 -O.OOlOi 0.9233 0.3821 
V 0.0223 - 0.0030i -0.3814 0.9241 / 



(5a) 



(56) 



where we have simply assumed a value for the electron-neutrino mass to which our analysis 
is not very sensitive and constructed the lepton mixing matrix by making use of the unitarity 
conditions with the same phase in (5b) as in (3b). 

We evolve the masses and mixing matrices to the SUSY GUT scale and use the extended 
Kusenko[|J procedure to construct the mass matrices numerically. The simplest SO(10) 
structure for the mass matrices is found with x q = and xe = 0.88, in which case the 
matrices have the following SO(10) transformation properties: 



M 



M N mrac _ dmc/(126; 126; 10) 
/ 10', 126 10', 126' 10' \ 



M 



D 



M 



10', 126' 



10' 



(6a) 
(66) 



126 

V 10' 10' 10 / 

Note that the same 10 is assumed to contribute to the 33 elements of the above mass 

matrices, with Yukawa coupling unification achieved for tan/? ~ 48.9. 

In this scenario we are able to construct a simple SO(10) model with just nine independent 
parameters for the following four matrices, such that 



M 



u 



diag(F', E', C) 



M N Dirac = diag{-?>F\ -3E', C) 





( ° 


A 






(F 





D\ 


M D = 


A 


E 


B 


M E = 





-3E 


B 




\D 


B 


cj 




U 


B 


c) 



(7a) 



where only D is complex and 

C'/C = v u /v d , 



E'/E = -AF'/F = w u /w d 



(76) 



in terms of the ratios of the 10 and of the 126 vacuum expectation values associated with 
the diagonal Yukawa couplings. In this model for scenario (A), two 10's and two 126's are 
required as indicated in (6a,b), while four texture zeros appear in the two matrices for M u 
and M D and for M N ° irac and M E . 



With 



F' = -1.098 x 10" 3 , 

C = 2.4607, 

E = -0.3830 x 10 _1 , 

F = -0.5357 x 10- 3 , 

A = -0.9700 x 10- 2 , 



E' = 0.314, C = 120.3 (8a) 

so v u /vd = tan (3 = 48.9 
hence w u /wd = —8.20 

(86) 

B = 0.8500 x lO" 1 
D = (0.4200 + 0.4285i) x 10~ 2 
in GeV, the masses and mixing matrices are calculated at the GUT scale by use of the 
projection operator technique of Jarlskog|| and then evolved to the low scales. The following 
low-scale results emerge for the quarks: 



m u (lGeV) =5.10MeV, 
m c (m c ) = 1.27 GeV, 
m t (m t ) = 150 GeV, 



/ 



V ( 



CRM 



0.9753 
-0.2209 



0.2210 
0.9747 



m d (lGeV) = 9.33 MeV 
m s (lGeV) = 181 MeV 
m b (m b ) = 4.09 GeV 

(0.2089 - 0.2242i) x 10- 2 \ 



(9a) 



(96) 



/ 



0.0444 

V 0.0078 - 0.0022i -0.0438 - 0.0005i 0.9994 
which are to be compared with the input starting data given in (3a) and (3b). 

In the absence of any VEVs coupling the left-handed neutrino fields together, we observe 

that the heavy righthanded Majorana neutrino mass matrix can be computed at the GUT 

scale from the seesaw mass formula 



which can be well approximated by the nearly geometric form 



(10a) 



M 



R 



F" 

-IVf^W' 



- 1 ^F JF W i - 1 \/F ir C"e i *v" \ 



E" 



\-\ \/F"C"e^v" - 1 VE"C"e^B" 



C" J 



(106) 
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where E" = \^F"C" and <f)B" = — <Pd"/3. In terms of the three additional parameters 
C" = 0.6077 x 10 15 , F" = 0.1745 x 10 10 and D » = 45°, the resulting heavy Majorana 
neutrino masses are determined to be 



(10c) 



M Rl = 0.249 x 10 9 GeV 
M R2 = 0.451 x 10 12 GeV 
M Ra = 0.608 x 10 15 GeV 

From the model parameters, the seesaw formula and Jarlskog's projection operator 
technique, || the light lepton masses and their mixing matrix can be constructed at the 
GUT scale and then evolved downward to the low scales where we find 



m 
m 
rn 



0.534 x 10" 5 eV, m 
0.181 x 10~ 2 eV, m 
0.135 eV, m T 



0.504 MeV 

105.2 MeV (Ha) 
1.777 GeV 



and 



V, 



LEPT 



( 0.9990 0.0451 
-0.0422 0.9361 
V 0.0174 - 0.0024i -0.3799 - O.OOOli 



-0.029 - 0.227i) x 10~ 2 \ 
0.3803 

0.9371 / 



(116) 



The agreement with the initial input values in (5a,b) is excellent. 



2.2 Lepton Scenario (B) involving the Non-adiabatic MSW Solar Depletion 
Effect and the Cocktail Model for Mixed Dark Matter 



In this scenario, the parameters of the 23 mixing plane given in (4) are replaced by 



5m 2 23 = 49 eV 2 , sin 2 2 8 23 = 10" 



(12a) 



with the tau-neutrino now assumed to account for the 30% hot dark matter component of 
the cocktail model H for mixed dark matter with a mass of tji Ut = 7.0 eV. The mixing angle 



V, 



(B) 



LEPT 



has been set close to the present upper bound from accelerator experiments, so the mixing 
matrix is now given by 

/ 0.9990 0.0447 (-0.289 - 0.129i) x 10~ 2 \ 

-0.0446 0.9989 0.0158 (126) 

\0.0036- 0.0013i -0.0157 0.9998 J 

Following the same general procedure as applied for scenario (A), we find the simplest 
SO (10) structure is obtained for x q = 0.5 and xe = 0.0, but only one texture zero is then 
present. Instead we make use of the bases where x q = and xg = 0.3 as in scenario (A). 
The quark mass matrices are then exactly the same as before, while the SO (10) structures 
are more complicated as seen from 



M u ~ M Noirac ~ diag(10, 126; 126; 10) 



/ 10, 10', 126 10', 126' 10' \ 



M 



D 



M 



E 



(13a) 



(136) 



10', 126' 126 10', 126' 

V 10' 10', 126' 10 J 
Without repeating all the details, we simply write down the form of the model matrices 

derived in this scenario and find 



M u = diag(F\ E' , C) 
/0 A D\ 

M D = A E B 

K D B C ) 



M Nl 



diag(-2.5F' , -3E', C) 



{ -F 

1 6 



M l 



-\A D\ 



-\A -3E IB 
\ D \B C ) 



(14) 



where again only D is complex. In this model for scenario (B), two 10's and two 126's 
are required as indicated in (13a,b), while four texture zeros appear in the two matrices for 
M u and M D but only three texture zeros appear for M Nmrac and M E . Since the quark 
mass matrices must lead to the same numerical results as in scenario (A), the values for the 
parameters introduced above are those given in (8a) and (8b). 
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The heavy right-handed Majorana mass matrix can be approximated by the nearly geo- 
metric form 



M R = 



F" 



-2VF"E" -\/F ir C r 'e i ^" \ 
2\/T jr E" -E" l^fWcy'e^B" 

y-y/F^e^' 2y/E 7i C JJ e i ' t> B" C" ) 



(15a) 



where E" = l^F^C", aside from an overall sign. With C" = 0.2323 x 10 14 , F" = 0.1096 x 
10 11 and 4>D" — 18.7°, 4>b" = 23.0° and <p C " = 41.8°, we find the resulting heavy Majorana 
neutrino masses for this case are determined to be 

M Rl = 0.841 x 10 9 GeV 

M R2 = 0.312 x 10 12 GeV (156) 
M Ra = 0.235 x 10 14 GeV 

By again making use of the simplified matrices at the GUT scale first to compute the 
lepton masses and mixing matrix Vlept by the projection operator technique of Jarlskog 
and then to evolve the results to the low scales, we find at the low scales for the (B) scenario 



m Ue = 0.544 x 10" 6 eV, 
m„„ = 0.242 x 10" 2 eV, 
m UT = 6.99 eV, 



m 
m 



i' 



0.511 MeV 
107.9 MeV 



(16a) 



m T = 1.776 GeV 



and 



V, 



LEPT 



( 0.9992 0.0410 (0.150 - 0.107i) x 10~ 2 \ 

-0.0411 0.9991 0.0113 



\ -0.0010 -O.OOlli -0.0123 0.9999 
to be compared with the initial low scale input in (5a) and (12). 



(166) 



3 Summary 

In this talk we have sketched a procedure which enables one to construct fermion mass 
matrices at the GUT scale which yield the low energy data taken as input. The models 
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constructed for the two neutrino scenarios work well in the SO(10) SUSY GUT framework 
with relatively few parameters, but the structure exhibited for scenario (B) with a 7 eV 
tau-neutrino is not as simple as that for scenario (A) based on the observed muon-neutrino 
atmospheric depletion effect. Discrete family symmetries giving rise to these models are now 
under investigation. Our general procedure to construct mass matrices can be applied to 
other symmetry-based frameworks as well. 

The research of CHA was supported in part by Grant No. PHY-9207696 from the Na- 
tional Science Foundation, while that of SN was supported in part by the U.S. Department 
of Energy, Grant No. DE-FG05-85ER 40215. 
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